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\ We report the results of calculations of pseudoscalar-isoscalar hadronic current correlators us- 
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ing the Nambu-Jona-Lasinio model and the real-time finite-temperature formalism. (This work 



■ represents a continuation of results reported previously for other current correlators.) Results 

> . 

£■ — , are presented for the temperatures range 1.2 < T/T c < 6.0, where T c is the temperature of the 

00 

confinement-deconfinement transition, which we take to be T c = 170 MeV. Some resonant features 



are seen in our calculations. In order to understand the origin of these resonances, we have per- 
formed relativistic random phase approximation (RPA) calculations of the temperature-dependent 
spectrum of the r\ mesons for T < T c . For the RPA calculations, use is made of a simple model 



in which we introduce temperature-dependent constituent quark masses calculated in a mean-field 
approximation and a temperature-dependent confining interaction whose form is motivated by re- 



cent studies made using lattice simulations of QCD with dynamical quarks. We also introduce 
temperature-dependent coupling constants in our generalized NJL model. Our motivation in the 
latter case is the simulation of the approach to a weakly interacting system at high temperatures 
and the avoidance of ll n condensates" which would indicate instability of the ground state of the 
model. We present some evidence that supports our use of temperature-dependent coupling con- 
stants for the NJL model. We suggest that our results may be of interest to researchers who use 
lattice simulations of QCD to obtain temperature dependent spectral functions for various hadronic 
current correlation functions. 
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I. INTRODUCTION 



Recently we have seen a good deal of interest in the calculation of hadronic current corre- 
lation functions at finite temperature using lattice simulations of QCD [1-5]. The procedure 
involves the calculation of the correlation function in Euclidean space and the application 
of a generalized Laplace transform to obtain the spectral functions. Various temperature- 
dependent resonant structures are found in the spectral functions. It is of interest to obtain 
further understanding of the nature of these resonances using chiral Lagrangian models, 
such as that of the Nambu-Jona-Lasinio (N JL) model [6] , that are successful in reproducing 
some of the low-energy properties of QCD. We have made some calculations of spectral 
functions making use of a generalized Nambu-Jona-Lasinio model [7]. In that work we 
studied correlators for currents with the quantum numbers of the ao, fo, 7r, and p mesons. 
In the present work we extend our study to the pseudoscalar-isoscalar correlators. In this 
case the correlators are calculated in terms of excitations that have the quantum numbers 
of the r] mesons. One of our goals in the present study is to relate the pseudoscalar-isoscalar 
spectral functions to the properties of the 77 mesons at finite temperature. The temperature 
dependence of the spectrum of the eta mesons is obtained using our generalized NJL model 
which includes a covariant model of confinement [8-13]. At zero temperature our model 
provides an excellent fit to the properties of the 7/(547), t/(958) mesons and their radial ex- 
citations [12]. (The model gives quite satisfactory results for the decay constants and mixing 
angles of the 77 and 7/ mesons.) In order to obtain the spectrum of the eta mesons at finite 
temperature, we make some assumptions which lead to relatively simple calculations. For 
example, we use temperature-dependent mass values for the constituent quarks obtained in 
a mean-field analysis and also introduce a temperature-dependent confining interaction that 
is motivated by some recent results obtained in lattice simulations of QCD with dynamical 
quarks. In addition, we have introduced a temperature dependence of the coupling constants 
of the NJL model. Our original motivation for introducing that temperature dependence 
was to model the various physical mechanisms that work against the development of pion 
condensation [14]. However, in the present work we argue that, if the NJL model is to be 
used at high temperatures, such temperature dependence is needed to avoid results that 
are inconsistent with what is known about QCD thermodynamics. For our studies we have 
introduced coupling constants G(T) = G(0)[1 — 0.17 T/T c ], where T c is the temperature 
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of the confinement-deconfinement transition, which we take to be T c = 170 MeV. At that 
temperature the coupling constants are reduced by 17%. The coupling constants are equal 
to zero beyond T = 5.88 T c . Thus, we see that in this model the short-range interaction 
is present in the range < T < 5.88 T c , while the confining interaction is taken to vanish 
for T > 1.2 T c . That feature of the model is also consistent with what is known concerning 
QCD thermodynamics. 

Ideally, it would be preferable if we could calculate the temperature-dependent spectrum 
of the rj mesons, including the radial excitations, in the imaginary-time Matsubara formalism 
or in the real-time finite-temperature formalism [15, 16]. That is a formidable task which is 
beyond the scope of the present work, but might be considered at some future time. (On the 
whole, there are some unresolved questions concerning the use of the field-theoretic finite- 
temperature theory in the confined phase of QCD which are not present when considering 
the deconfined phase.) In those cases in which our results may be compared to those of 
the Matsubara formalism we find general agreement [17]. However, we need to assume 
for the present work that, with the use of temperature-dependent constituent quark mass 
values, a temperature-dependent confining interaction and temperature-dependent coupling 
constants, we have introduced the most important features that influence the 77 spectrum 
for T < T c . 

The organization of our work is as follows. In Section II we discuss the calculation of 
polarization functions of the NJL model at finite temperature. In Section III we describe the 
calculation of the pseudoscalar-isoscalar hadronic current correlation functions and present 
various results of our numerical computations. Since the use of temperature-dependent cou- 
pling constants is an unusual feature of our model, we provide some justification for the use 
of such coupling constants in Section IV. In Section V we introduce the Lagrangian of our 
model that is used for our RPA calculations. We make reference to the RPA equations that 
were used to calculate the properties of the 77 mesons at T = in earlier work and describe 
the motivation for the introduction of a temperature-dependent confining interaction. We 
present some results of our numerical calculations and relate the results of our RPA calcu- 
lations for T < T c to the results obtained for the correlation functions for T > T c . Finally, 
Section VI contains some further discussion and conclusions. 
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II. POLARIZATION FUNCTIONS AT FINITE TEMPERATURE 

In an earlier work we carried out a Euclidean-space calculation of the up, down, and 
strange constituent quark masses taking into account the 't Hooft interaction and our con- 
fining interaction [18]. The 't Hooft interaction plays only a minor role, but does provide 
coupling of the equations for the various constituent masses. If we neglect the confining in- 
teraction and the 't Hooft interaction in the mean field calculation of the constituent masses, 
we can compensate for their absence by making a modest change in the value of Gs, the 
coupling constant of the NJL model. For the calculations of this work we calculate the quark 
masses using the formalism presented in the Klevansky review [19]. (Note that our value of 
G s is twice the value of G used in that review.) The relevant equation is Eq. (5.38) of Ref. 
[19]. Here, we put \i = and write 



where A = 0.631 GeV is a cutoff for the momentum integral, (3 = l/T and E p — [p 1 + 
m 2 (T)]V 2 . i n our calculations we replace G by G s (T)/2 and solve the equation 



with G S (T) = 11.38[1 - 0.17T/TJ GeV, m° u = 0.0055 GeV and m° s = 0.130 GeV. Thus, we 
see that Gs(T) is reduced from the value Gs(0) by 17% when T = T c . The results obtained 
in this manner for m u (T) and m s (T) are shown in Fig. 1. Here, the temperature dependence 
we have introduced for Gs{T) serves to provide a somewhat more rapid restoration of chiral 
symmetry than that which is found for a constant value of Gs- That feature and the 
temperature dependence of the confining potential leads to the deconfinement of the light 
mesons considered here at T > T c . 

For the calculation of polarization functions for T > 1.2 T c we may neglect the confining 
interaction. However, we include temperature-dependent quark mass values in our calcula- 
tions. The basic polarization functions that are calculated in the NJL model are shown in 
Fig. 2. We will consider calculations of such functions in the frame where P = 0. In our 
earlier work, calculations were made after a confinement vertex was included. That vertex is 
represented by the filled triangular region in Fig. 2. However, we here consider calculations 
for T > 1.2T C where confinement may be neglected. 




(2.1) 




(2.2) 
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FIG. 1: Temperature dependent constituent mass values, m u (T) and m s (T), calculated using Eq. 
(2.2) are shown. Here m° u = 0.0055 GeV, m° s = 0.130 GeV, and G{T) = 5.691 [1 - 0.17(T/T C )], if 
we use Klevansky's notation [19]. 




FIG. 2: The upper figure represents the basic polarization diagram of the NJL model in which 
the lines represent a constituent quark and a constituent antiquark. The lower figure shows a 
confinement vertex [filled triangular region] used in our earlier work. For the present work we 
neglect confinement for T > 1.2T C , with T c = 170 MeV. 

The procedure we adopt is based upon the real-time finite-temperature formalism, in 
which the imaginary part of the polarization function may be calculated. Then, the real 
part of the function is obtained using a dispersion relation. The result we need for this 
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work has been already given in the work of Kobes and Semenoff [20]. (In Ref. [20] the 
quark momentum in Fig. 2 is k and the antiquark momentum is k — P. We will adopt that 
notation in this section for ease of reference to the results presented in Ref. [20].) With 
reference to Eq. (5.4) of Ref. [20], we write the imaginary part of the scalar polarization 
function as 

1 - JS ^ = J^e^ (2.3) 

{(1 - ni (k) - n 2 (k))5(p° - E 1 (k) - E 2 {k)) 
-(ni(fc) - n 2 (k))5(p° + E^k) - E 2 {k)) 

-M*) - ni(fc))<f(P° - E^k) + E 2 {k)) 
_(1 _ ni (k) - n 2 (k))5(p° + E 1 (k) + E 2 (k))} . 

Here, E^k) = [k 2 + m^T)] 1 / 2 . Relative to Eq. (5.4) of Ref. [20], we have changed 
the sign, removed a factor of g 2 and have included a statistical factor of 2N C , where the 
factor of 2 arises from the flavor trace. In addition, we have included a Gaussian regulator, 
exp[— k 2 1 a 2 ], with a = 0.605 GeV, which is the same as that used in most of our applications 
of the NJL model in the calculation of meson properties [8-15]. We also note that 

and 

n2 ^ = e PE 2 l) + x • ( 2 - 5 ) 

For the calculation of the imaginary part of the polarization function, we may put k 2 = 
m\{T) and (k — P) 2 = m 2 {T), since in that calculation the quark and antiquark are on- 
mass-shell. We will first remark upon the calculation of scalar correlators. In that case, the 
factor p s m Eq. (2.3) arises from a trace involving Dirac matrices, such that 

Ps = -TVP + mi )(ft-p + m 2 )\ (2.6) 
= 2P 2 -2(m! + m 2 ) 2 , (2.7) 

where mi and m 2 depend upon temperature. In the frame where P = 0, and in the case 
mi = m 2) we have Ps — 2-P 2 (l — 4m 2 /P 2 ). For the scalar case, with m x = m 2) we find 

Im J S (P 2 ,T) = ^ (l " ^) 3/2e " P/a2 [ 1 - 2»i(*)l , (2-8) 
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where 

k 2 = ^j--m 2 (T). (2.9) 

We may evaluate Eq. (2.8) for m(T) = m u (T) = rrid(T) and define Im J^(P 2 ,T). Then 
we put m(T) = m s (T), we define Im Jf (P 2 , T). These two functions are needed for a 
calculation of the scalar-isoscalar correlator. (Note that the factor of 2 arising from the 
flavor trace should be moved when we define the polarization function for a specific quark 
flavor.) The real parts of the functions Jf (P 2 ,T) and jf (P 2 ,T) may be obtained using a 
dispersion relation, as noted earlier. 

For pseudoscalar mesons, we replace [3s by 

P P = -Ti[i l5 {lt + mi)i 75 (^ - P + m 2 )\ (2.10) 
= 2P 2 - 2(mi - m 2 ) 2 , (2.11) 

which for mi = m 2 is f3p = 2P 2 in the frame where P = 0. We find, for the tt mesons, 

toJ'(P>,T) = ^ (l- ^) 1/2 e- t2 ^[l-2n l( m, (2.12) 

where k 2 = Pq/4 — m 2 (T), as above. Thus, we see that, relative to the scalar case, the 
phase space factor has an exponent of 1/2 corresponding to a s- wave amplitude, rather than 
the p-wave amplitude of scalar mesons. For the scalars, the exponent of the phase-space 
factor is 3/2, as seen in Eq. (2.8). 

For a study of vector mesons we consider 

/%, = Tr[ 7M (^ + - P + m 2 )\ , (2.13) 

and calculate 

gTPlu = 4[P 2 - ml - mj + 4 mi m 2 ] , (2.14) 

which, in the equal- mass case, is equal to 4Pq + 8m 2 (T), when mi = m 2 and P = 0. Note 
that for the elevated temperatures considered in this work m u {T) = m^(T) is quite small, 
so that 4P 2 + 8m 2 (T) can be approximated by 4P 2 when we consider the p meson. The 
generalization of these results for the study of the pseudoscalar-isoscalar correlators will be 
taken up in the next section. 
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III. CALCULATION OF HADRONIC CURRENT CORRELATION FUNCTIONS 

In this section we consider the calculation of temperature-dependent hadronic current 
correlation functions. The general form of the correlator is a transform of a time-ordered 
product of currents, 



where the double bracket is a reminder that we are considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp,i{x) = 
q(x)i'j5\ l q(x), where, in the case of the n mesons, i = 1,2, and 3. For the study of 
pseudoscalar-isoscalar mesons, we again introduce jp,i(x) = q(x)X l q(x), but here i = 
for the flavor-singlet current and i = 8 for the flavor-octet current. 

In the case of the n mesons, the correlator may be expressed in terms of the basic vacuum 
polarization function of the NJL model, J P (P 2 ,T) [19, 21, 22]. Thus, 



where G n (T) is the coupling constant appropriate for our study of the tt mesons. We have 
found G n (0) = 13.49 GeV" 2 by fitting the pion mass in a calculation made at T = 0, with 



The calculation of the correlator for pseudoscalar-isoscalar states is more complex, since 
there are both flavor-singlet and flavor-octet states to consider. We may define polarization 
functions for u, d and s quarks: J U (P 2 ,T), J d (P 2 ,T) and J S (P 2 ,T). (We recall that the 
factor of 2 arising from the flavor trace is not included when these functions are calculated.) 
In terms of these polarization functions we may then define 




(3.1) 




(3.2) 



m u = md = 0.364 GeV [14]. 




(3.3) 




(3.4) 



and 



±[J u (P 2 ,T) + J d (P 2 ,T)+4J s (P 2 ,T)]. 



(3.5) 
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We also introduce the matrices 



J(P 2 ,T) 



G{T) = 



J m (P\T) J s(P 2 ,T) 
J 80 (P 2 ,T) J 88 (P 2 ,T) 



Goo(T) G 08 (T) 
G 80 (T) G 88 (T) 



(3.6) 



(3.7) 



and 



C(P 2 ,T) 



(3.8) 



(3.9) 



(3.10) 



C o(P 2 ,T) C 08 (P 2 ,T) 
^C 80 (P 2 ,T) C 88 (P 2 ,T) 

We then write the matrix relation 

C(P 2 ,T) = J(P 2 ,T)[1 - G{T)J(P 2 ,T)\- 1 

For some purposes it may be useful to also define a t matrix 

t(P 2 ,T) = [1 - G{T)J{P 2 ,T)]- 1 G{T) , 

where t(P 2 ,T) has the structure shown in Eqs. (3.6)-(3.8). The same resonant structures 
are seen in both C{P 2 ,T) and t(P 2 ,T). 

Some of our results for the imaginary parts of the pseudoscalar-isoscalar correlators 
Coo(P 2 ), C 88 (P 2 ) and C 08 (P 2 ) are shown in Figs. 3, 4 and 5, respectively. In these fig- 
ures the values are T/T c = 1.2 [solid line], T/T c = 1.6 [dashed line], T/T c = 2.0 [dotted 
line], T/T c = 4.0 [dashed-dotted line] and T/T c = 6.0 [dashed- (double) dotted line]. There 
is a large peak seen in Figs. 3-5 at about 775 MeV. It is worth noting that the state that 
evolves from the ?/(958) with increasing temperature has a mass of about 750 MeV for 
T ~ T c . However, an analysis of the mixing angle for the state at 775 MeV shows that 
it is mainly an ss state. Further work is needed to understand the relation between the 
bound-state spectrum for T < T c and the resonant structures seen for T > T c . We discuss 
the temperature dependence of the r\ spectrum in the Section V. 



IV. TEMPERATURE-DEPENDENT COUPLING CONSTANTS OF THE NJL 
MODEL 

Since the introduction of temperature-dependent coupling constants for the NJL model is 
a novel feature of our work, we provide arguments in this section to justify their introduction. 
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FIG. 3: The imaginary part of the pseudoscalar-isoscalar correlator Cqq{P 2 ) is shown. Here, 
T/T c = 1.2 [solid line], 1.6 [dashed line], 2.0 [dotted line], 4.0 [dashed-dotted line] and 6.0 [dashed- 
(double) dotted line]. In this work we use G 00 = 8.09 GeV~ 2 , G 88 = 13.02 GeV -2 and G 08 = 
-0.4953 GeV" 2 . 
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FIG. 4: The imaginary part of the correlator Cgs{P 2 ) is shown. [See caption to Fig. 3.] 
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FIG. 5: The imaginary part of the correlator Co 8 (-P 2 ) is shown. [See caption to Fig. 3.] 

We make reference to Fig. 1.3 of Ref. [16]. That figure shows the behavior of the ratio e/T 4 
and 3P/T 4 for the pure gauge sector of QCD. Here e is the energy density and P is the 
pressure. Ideal gas behavior implies e = 3P. The values of e/T 4 and 3P/T 4 are compared 
to the value 6sb /T 4 = 87r 2 /15 for an ideal gluon gas. It may be seen from the figure that 
at T = 3T C there are still significant differences from the ideal gluon gas result. Deviations 
from ideal gas behavior become progressively smaller with increasing T/T c and could be 
considered to be relatively unimportant for T/T c > 5. 

The use of our energy-dependent coupling constants is meant to be consistent with the 
approach to asymptotic freedom at high temperature. In order to understand this feature 
in our model, we can calculate the correlator with constant values of Goo, G 88 and Go 8 and 
with G 00 (T) = Goo[l - 0.17 T/TJ, etc. (In this work we use G 00 = 8.09 GeV~ 2 , G 88 = 13.02 
GeV- 2 and G 08 = -0.4953 GeV- 2 .) 

We now consider the values of Im G 88 (P 2 ) for T/T c = 4.0. In Fig. 6 we show the values 
of Im G 88 (P 2 ) calculated in our model with temperature-dependent coupling constants as a 
dashed line. The dotted line shows the values of the correlator for Goo — G 88 = Go 8 = 0, 
while the solid line shows the values when the coupling constants are kept at their values 
at T = 0. We see that we have some resonant behavior in the case the constants are 
temperature independent. 
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FIG. 6: The imaginary part of the correlator CssiP 2 ) is shown for T/T c = 4.0. The dashed line 
is the result for the temperature-dependent coupling constants of our model, while the solid line 
represents the results for coupling constants kept at their T = values. [See caption to Fig. 3.] 
The dotted line shows the values of the correlator when the coupling constants are set equal to 
zero. 

In Fig. 7 we show similar results for T/T c = 5.88. Here the temperature-dependent 
coupling constants are equal to zero, so that the lines corresponding to the dashed and 
dotted lines of Fig. 6 coincide. The solid line again shows some resonant behavior at a 
value of T/T c , where we expect only very weak interactions associated with asymptotic 
freedom. We conclude that the model with constant values of the coupling constants yields 
unacceptable results, while our model, which has temperature-dependent coupling constants, 
behaves as one may expect, when the results of lattice simulations of QCD thermodynamics 
are taken into account. 
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FIG. 7: The imaginary part of the correlator Css{P 2 ) is shown for T/T c = 5.88. [See caption to 
Fig. 6.] Here the dashed and dotted lines of Fig. 6 coincide. 

V. CALCULATION OF MESON PROPERTIES AT FINITE TEMPERATURE IN 
A GENERALIZED NJL MODEL WITH CONFINEMENT 

It is useful to record the Lagrangian used in our calculations of meson properties 
C = q(i<? - m°)q + ^ ^[(gA'g) 2 + (qi^Xq) 2 } 

i=0 

-^E[(^7^) 2 + (^757,g) 2 ] 

i=0 

+ ^{det[g(l + 75 )g] + det[g(l - l5 )q}} 

(5.1) 

Here, m° is a current quark mass matrix, m° = diag(m°, mjj, m° s ). The X 1 are the Gell- 
Mann (flavor) matrices, A = a/2/3 1 with 1 being the unit matrix. The fourth term on the 
right-hand side of Eq. (5.1) is the 't Hooft interaction. Finally, C con f represents the model 
of confinement we have used in our work. 

As noted earlier, we have recently reported results of our calculations of the temperature 
dependence of the spectra of various mesons [14]. These calculations were made using our 
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generalized NJL model which includes a covariant model of confinement. We have presented 
results for the 7r, K, a , fo and Kq mesons in Ref. [14]. The equations that we solve are 
of the form of relativistic random-phase-approximation equations. The derivation of these 
equations for pseudoscalar mesons is given in Ref. [13], where we discuss the equations for 
pionic, kaonic and eta mesons. The equations for the eta mesons are the most complicated, 
since we consider singlet-octet mixing as well as pseudoscalar-axial-vector mixing. In that 
case there are eight vertex functions to consider, rp^, T\^, Tpg, Tj^g, Tp^, F^q, r^g, 
where P refers to the 75 vertex and A refers to the 7075 vertex, which mixes with the 
75 vertex. Corresponding to the eight vertex functions one may define eight wave function 
amplitudes [13]. Since the RPA equations for the study of the eta meson are quite lengthy 
[13], we do not reproduce them here. 

In the RPA equations we replace m u and m s by m u (T) and m s {T) of Fig. 1 and use 
the temperature-dependent coupling constants described earlier in this work. In addition 
to the temperature dependence of the coupling constants and constituent mass values, we 
also introduced a temperature-dependent confining potential, whose form was motivated by 
recent lattice simulations of QCD in which the temperature dependence of the confining 
interaction was calculated with dynamical quarks [23]. (See Fig. 8.) In order to include 
such effects, we modified the form of our confining interaction, V c (r) = M~exp[— fj,r], by 
replacing \i by 

m = 7 — , ,n , (5-2) 



rj-\ \ - 



with /jl = 0.010 GeV. The maximum value of V c (r,T) is then 

^(T) = -^-, (5.3) 

= «[l-0.7(T/r e ) 2 ] 
A*oe 

with r max = l//i(T). To better represent the qualitative features of the results shown in Fig. 
8, we use V c (r,T) = kt exp[-/x(T)r] for r < r max and V c (r,T) = V£ ax (T) for r > r max . 
We also note that we use Lorentz-vector confinement and carry out all our calculations in 
momentum-space. (The value of n used in our work is 0.055 GeV 2 .) Values of V c (r, T) are 
shown is Fig. 9. 
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FIG. 8: A comparison of quenched (open symbols) and unquenched (filled symbols) results for the 
interquark potential at finite temperature [23]. The dotted line is the zero temperature quenched 
potential. Here, the symbols for T = 0.80T C [open triangle], T = 0.88T C [open circle], and T = 
0.94T C [open square], represent the quenched results. The results with dynamical fermions are given 
at T = 0.68T C [solid downward-pointing triangle], T = 0.80T C [solid upward-pointing triangle], 
T = 0.88T C [solid circle], and T = 0.94T C [solid square]. 

In Fig. 10 we show the results of our calculations of the temperature-dependent spectrum 
of the rj mesons. We show the behavior of the 77(547), 77'(958) and seven states which 
represent radial excitations. The energies of the additional states found when diagonalizing 
the RPA Hamiltonian are represented by dots in the range < T/T c < 0.6. We note that 
the masses of the nodeless states (the r/ and rj') are fairly constant over a broad range of 
temperatures. That characteristic seems to be a feature of the behavior of pseudo-Goldstone 
bosons at finite temperature. 

As the temperature is increased, fewer states are bound by the confining field which 
decreases in magnitude with increasing temperature. At T = T c only the state that evolves 
from the r/(547) is bound. That state disappears from the spectrum for T > T c . We believe 
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FIG. 9: The potential V c (r,T) is shown for T/T c = [solid line], T/T c = 0.4 [dotted line], T/T c = 
0.6 [dashed line], T/T c = 0.8 [dashed-dotted line], T/T c = 0.9 [short dashes], T/T c = 1.0 [dashed- 
(double) dotted line]. Here, V c (r,T) = nrexp[-fji(T)r], with fi(T) = 0.01GeV/[l - 0.7(T/T C ) 2 ] 
and k = 0.055 GeV 2 . 

that the crossing of levels seen at T/T c = 0.9 is due to the rapid decrease of m u (T) relative 
to m s (T) with increasing T/T c . [See Fig. 1.] That feature could lead to a (predominately) 
nn state with a node to have a lower energy than a (predominantly) nodeless ss state. 

VI. DISCUSSION AND CONCLUSIONS 

We believe is of interest to supplement lattice studies of hadronic current correlation 
functions with calculations made using chiral Lagrangian models of the type considered 
in this work. We have made some progress in exhibiting results for such correlators in 
Ref. [7] and in the present study. It might be of some interest to compare our results 
for our temperature-dependent RPA calculations with results obtained in the imaginary- 
time or real-time formalisms, if these formalisms could be modified so that calculations 
could be made in the confined phase of QCD. The study of radial excitations in these finite- 
temperature theories may be quite difficult since their study requires a model of confinement. 
It is also very difficult to obtain information concerning radial excitations, if an analytic 
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T/T c 

FIG. 10: The temperature-dependent spectrum of the r\ mesons is shown. For the most highly 
excited radial excitations we represent the mass values obtained by dots. There are no bound 
states for T > T c . 

continuation to real time is necessary. 

One interesting feature of our analysis is the use of temperature-dependent coupling 
constants in the NJL model. In the present work, we have provided some justification for 
the introduction of such constants. Our work suggests that the coupling constants of the 
NJL model may also be density-dependent, since one expects that high density may play a 
similar role as high temperature, leading ultimately to a weakly interacting system at high 
density. We have introduced density-dependent coupling constants in Ref. [24] where we 
considered the confinement-deconfinement transition in the presence of matter. Since the 
study of matter at high density is a topic of active investigation [25-31], our suggestion of 
density-dependent coupling constants may have important consequences for such studies. 
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